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ABSTRACT 



Context. We present modelling of theySCep star HD 180642 based on its observational properties deduced from CoRoT and ground- 
based photometry as well as from time-resolved spectroscopy. 

Aims. We investigate whether present-day state-of-the-art models are able to explain the full seismic behaviour of this star, which has 
extended observational constraints for this type of pulsator. 

Methods. We constructed a dedicated database of stellar models and their oscillation modes tuned to fit the dominant radial mode 
frequency of HD 180642, by means of varying the hydrogen content, metallicity, mass, age, and core overshooting parameter. We 
compared the seismic properties of these models with those observed. 

Results. We find models that are able to explain the numerous observed oscillation properties of the star, for a narrow range in mass 
of 11.4-11.8 Mq and no or very mild overshooting (with up to 0.05 local pressure scale heights), except for an excitation problem of 
the £ = 3, pi mode. We deduce a rotation period of about 13 d, which is fully compatible with recent magnetic field measurements. 
The seismic models do not support the earlier claim of solar-like oscillations in the star. We instead ascribe the power excess at high 
frequency to non-linear resonant mode coupling between the high-amplitude radial fundamental mode and several of the low-order 
pressure modes. We report a discrepancy between the seismic and spectroscopic gravity at the 2.5o" level. 
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On the basis of observations by the small satellites WIRE 
(Buzasi 2004) and MOST (Matthews 2007), high-precision 
time-resolved space photometry has been obtained by the ongo- 
ing missions CoRoT (Baglin et al. 2006; Auvergne et al. 2009) 
and Kepler (Burocki et al. 2010; Gilliland et al. 2010) provid- 
ing significant progress in particular subfields of asteroseismol- 
ogy. Notable discoveries have been made of oscillation proper- 
ties unknown prior to the era of space asteroseismology for solar- 
like oscillations, not only in both FGK-type main-sequence stars 
(e.g., Michel et al. 2008, Chaplin et al. 201 1) and red giants (e.g., 
De Ridder et al. 2009, Beck et al. 201 1, Bedding et al. 201 1), but 
also in a massive O-type star (Degroote et al. 2010b). 

Thanks to the continuity between space- and ground-based 
observations, the discovery of hybrid pulsators in various evolu- 
tionary stages, i.e., stars pulsating in pressure (p) and gravity (g) 
modes simultaneously, will also shed new light on the details of 
stellar structure models. While a few of these hybrid pulsators 
had been known prior to the space missions, large amounts of 
them were found from the space photometry, all along the main 
sequence from spectral type B (e.g., Balona et al. 2011) to F 
(Grigahcene et al. 2010). Seismic diagnostics suitable for inter- 
preting these cases in terms of modes probing both the core re- 
gion and the outer envelope are yet to be developed, but it is 
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clear that the numerous gravity modes offer a unique opportu- 
nity to test the physics in the deep interior regions adjacent to 
the core. Initial efforts to explore this region concluded that the 
input physics of the current models has to be improved to accu- 
rately describe chemically inhomogeneous regions in order to be 
compliant with the space data (Degroote et al. 2010a). 

For the most massive pulsators of spectral types O and BO-2, 
space photometry so far provided more questions than answers. 
While the detection of pulsations in photometric data was finally 
established for O-type stars, there seems to be a different cause 
of the pulsations. The 09V star HD 46202 behaves in a similar 
way to known ySCep stars (e.g., Aerts et al. 2010 for an exten- 
sive overview of all types of pulsators across the Hertzsprung- 
Russell diagram) but none of its modes is predicted to be excited 
by the theory that is able to explain the ySCep instability strip 
(Briquet et al. 2011). Moreover, the 08.5V star HD 46149 ap- 
pears to oscillate in stochastic modes fulfilling scaling laws and 
echelle diagrams with clear ridges for modes of equal degree 
but consecutive radial orders as those observed for Sun-like stars 
(Degroote et al. 2010b). This implies that this star must have an 
outer convection zone of considerable size, a suggestion made 
for massive stars in general, independently of seismic data, from 
theoretical work (Cantiello et al. 2009). Three hotter O stars do 
not show any coherent classical pulsational signal, but rather an 
excess power at low frequency resembling red noise (Blomme et 
al. 2011). More accurate knowledge of the physics in the outer 
envelope of O stars is needed to understand the variability be- 
haviour detected from CoRoT (there are no O stars in the Kepler 
FoV). 
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Table 1. Frequencies and mode identification of the highest am- 
pHtude modes of HD 180642. 



Frequency 

(d-') 


MI method 


e 


m 


i 

n 


(km s-') 


5.48694 


ampl.ratios 












7.35867 


ampl.ratios 


Oor3 








8.40790 


LPVs 


3 


+2 


82 


38 




or 


3 


±1 


31 


78 




or 


1 


±1 


15 


119 




or 


2 


±1 


14 


104 




or 


2 


±1 


76 


26 




or 


2 


±2 


48 


52 




or 


3 


±3 


67 


26 



6.26527 


ii^iii^ ... ... ... ... 

none ... 


6.32482 


none ... 


7.10353 


none ... 


7.25476 


none ... 


8.77086 


none ... 


0.29917 


none ... 


0.89870 


none ... 



Notes. MI method can be any of: photometric amplitude ratios 
(ampl.rations), line-profile variations (LPVs), not available (none). 



Few BO-2 pulsators with or without emission lines have been 
observed by CoRoT, while most of the B/Be stars monitored by 
Kepler are cooler. A remarkable finding is the outburst signal 
detected in the B0.5IVe star HD 49330 (Huat et al. 2009), a phe- 
nomenon that has not been observed to occur in any of the other 
observed Be stars, which are of later spectral type. Theoretical 
work is needed to understand and interpret the observed fre- 
quency patterns in the Be pulsators (Neiner et al. 2009, Diago 
et al. 2009, Gutierrez-Soto et al. 2009). Progress is hampered 
by the mathematical complexity of describing the oscillations in 
these flattened stars, induced by their rapid rotation, often close 
to their critical velocity. 

To date CoRoT has observed two slow rotators without emis- 
sion lines situated in the ySCep instability strip. The star B0.5V 
star HD 51756 turned out to be a rotational modulation vari- 
able rather than a pulsator (Papics et al. 2011). In this work, 
we present seismic modelling of the only known ySCep class 
member monitored by CoRoT: the B1.5II-III star HD 180642 
(V1449Aql). 



2. Summary of observational constraints 

The ySCep class member HD 180642 has a visual magnitude of 
8.3 and a dominant radial mode with V amplitude 39mmag 
(Aerts 2000). The additional observational characteristics of 
HD 180642 we used as input for the seismic modelling were 
taken from Degroote et al. (2009, Paper I) and Briquet et al. 
(2009, Paper II). Paper I was devoted to the analysis of the 
CoRoT light curve of the star. Three versions of a model de- 
scription for the light curve were considered: one with classi- 
cal prewhitening producing 127 significant frequencies, one with 
harmonics and combination frequencies of 11 independent fre- 
quencies, and one where time-dependent amplitudes and phases 
were allowed. After careful statistical evaluation taking into ac- 
count penalties for the number of free parameters, the best fit 
to the CoRoT light curve was found to be the model based on 
33 frequencies among which 1 1 were independent ones (which 
are repeated in Table[TJ, 3 were harmonics of the dominant fre- 



quency and 19 were combination frequencies (listed in Table 2 
of Paper I and not repeated here). 

The list of nine highest-amplitude independent stable fre- 
quencies contains the one of the dominant radial mode known 
prior to the CoRoT launch (Aerts 2000), eight frequencies in 
the range [6.1, 8.5] d' , and two frequencies below Id"' , which 
could either be due to g modes or connected to the rotation of the 
star (or both) and agree to within a factor three. 

These ten frequencies have amplitudes roughly between 0.8 
and 6mmag (Paper I), i.e., far below the 39mmag of the dom- 
inant mode. The frequency spectrum was found to be atypical 
of a ySCep star, in the sense that ten sum and nine difference 
frequencies were also detected, in addition to three harmonics 
of the dominant frequency. Moreover, several of the combina- 
tion frequencies were found to have locked phases, a character- 
istic expected for modes undergoing non-linear mode interaction 
with the large-amplitude dominant oscillation. The 1 1 indepen- 
dent largest-amplitude frequencies found in the space photom- 
etry are listed in Table[T] along with their mode identification 
obtained by Aerts (2000) and in Papers I and II. In this work, we 
consider the nine highest frequencies in that table for the mod- 
elling procedure, i.e., the two frequencies below Id"' are not 
used. The total range of frequencies related to dominant inde- 
pendent modes and their combinations covers [0.3, 27.5] d"' . 

In addition to these 33 frequencies, variations with time- 
dependent behaviour were found in the residual light curve, 
within the frequency range of the harmonics and the combi- 
nation frequencies of the heat-driven modes (Belkacem et al. 
2009; Paper I). Belkacem et al. (2009) determined a frequency 
spacing of 1.17d"'for the range above 11.23 d"' and below 
26 d "' and interpreted it in terms of solar-like oscillations caused 
by stochastic forcing due to convective motions in the outer 
stellar layers. In Paper I, a value for candidate frequency spac- 
ings was recomputed over the range [4.32, 25.9] d"' for the three 
types of prewhitening procedures and led to 1.05, 1.11, and 
1.56 d "' , with an uncertainty of ~ 0.02 d"' (the correspond- 
ing four spacing values are 13.5, 12.1, 12.9, and 18.0/iHz). The 
physical interpretation of these spacings remains unclear, given 
that their value depends on the way the CoRoT light curve is 
prewhitened. We return to this problem in Sect. 13. 71 

In Paper II, we presented the results of an extensive ground- 
based multicolour photometric and high-resolution spectro- 
scopic campaign that we have organised. This campaign has led 
to the detection of three of the CoRoT frequencies in the pho- 
tometry and nine in the spectroscopy, among then the fourth 
harmonic of the dominant frequency, which is not found in the 
space photometry. Empirical mode identifation, based on photo- 
metric amplitude ratios and on line-profile variability, led to the 
results in Table[T] In particular, the frequency 7.35867 d"' was 
identified with either an ^' = or ^ = 3 mode from amplitude 
ratios, while the moment method was used to deduce informa- 
tion on the mode wavenumbers {(, m), the inclination angle, and 
the equatorial velocity for the frequency 8.40790d "' . The seven 
possibilities listed in Table[T] occurred with a much higher prob- 
ability than any of the other options for (d, m) and we allow any 
of those seven solutions in the modelling rather than taking only 
the best one. 

The average spectrum of HD 180642 was also analysed in 
Paper II to determine the fundamental parameters of the star, 
as well as its current surface abundance pattern. This resulted 
in Teff = 24500 ± 1 OOOK, log g = 3.45 + 0.15, and Z e 
[0.0083, 0.0142], where the interval for Z covers the dependence 
on the microturbulent velocity that could not be easily fixed and 
led to a systematic uncertainty. 
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We also recall from Paper II that the overall line broadening 
amounts to 44 km s ' . This is the combined effect of line broad- 
ening caused by the projected rotation velocity and the pulsa- 
tions. It is thus an upper limit to v sin /. The dominant mode has 
a radial-velocity amplitude of 39 km s"' (Paper II), which, taking 
into account limb-darkening effects, translates the overall broad- 
ening value into an upper limit of ~ 30 km i"' for v sin /. 

Finally, we point out that Hubrig et al. (201 1) obtained and 
analysed spectropolarimetric time series observations taken with 
the SOFIN echelle spectrograph at the 2.56m Nordic Optical 
Telescope, after the discovery of a magnetic field in the star 
(Hubrig et al. 2009). From their 13 data points, three domi- 
nant peaks emerged in the periodogram, with the highest peak 
at 0.072 d ' . This periodicity was interpreted in the framework 
of a rigid rotator model for which the period of the magnetic 
field variation corresponds to a stellar rotation period of 13.9d, 
assuming a centred dipole. 



3. Modelling procedure 

The frequencies above Id' listed in Table[T]are typical of the 
low-order p and g modes in ySCep stars (e.g., Aerts et al. 2010). 
The rotation period deduced from the magnetic field measure- 
ments combined with realistic values for the radius of a y6 Cep 
star lead to equatorial surface rotation velocities that are low rel- 
ative to the critical velocity. Ignoring the effects of rotation in 
computing evolutionary models and their oscillation frequencies 
for the low-order zonal modes is justified in this situation. In 
practise, this means that we ignore the terms proportional to the 
square of the rotation frequency in all the equations, including 
that of hydrostatic support. 

For the stellar modelling, we followed a procedure similar to 
the one adopted by Ausseloos et al. (2004) but for an updated and 
more extensive grid of models. The latter was computed with 
the Code Liegeois d'Evolution Stellaire (CLES, Scuflaire et al. 
2008a), using the input physics described in Briquet et al. (201 1) 
which is not repeated here. An extensive evolutionary model grid 
covering the entire core-hydrogen burning phase along the main 
sequence (MS), i.e., from the zero-age main-sequence (ZAMS) 
to the terminal-age main-sequence (TAMS), for masses M from 
7.6 M0 to 20.0 Mq in steps of 0. 1 M©, hydrogen mass fractions X 
from 0.68 to 0.74 in steps of 0.02, metallicities Z from 0.010 un- 
til 0.018 in steps of 0.002 and core overshooting parameters Oow 
between zero and 0.5 in steps of 0.05 (expressed in local pres- 
sure scale heights) were considered. The probability of catching 
the star beyond the TAMS, i.e., in the Hertzsprung gap, is so 
low that it is not necessary to consider hydrogen-shell burning 
models. For each of the MS models along the 27500 tracks, a 
theoretical oscillation frequency spectrum ignoring the effects of 
rotation, v„f , was computed in the adiabatic approximation with 
the code LOSC (Scuflaire et al. 2008b), covering the range of the 
low-order (typically up to « = 7) p and g modes of degree from 
^ = up to ^ = 4. 

Following Ausseloos et al. (2004), we did not restrict the 
model parameters a priori by imposing limits derived from 
classical spectroscopy, because discrepancies between the seis- 
mically and spectroscopically derived fundamental parameters 
have been reported for a few massive stars in the literature, in 
particular for the gravity (e.g.. Briquet et al. 2007 for the /J Cep 
star flOph and Briquet et al. 2011 for the 09V-type pulsator 
HD 46202). Hence, we instead used the spectroscopic informa- 
tion as a posteriori constraints, after having fitted the nine high- 
est frequencies listed in Table[T] 



In the following, we considered models able to fit the domi- 
nant radial mode frequency, after which we continued with those 
for which this mode is excited. For the surviving models, we then 
required that the additional frequencies higher than the dom- 
inant mode frequency be fitted as well. Next, we considered 
the models fulfilling all these requirements, as well as the con- 
straint on the observed effective temperature. Finally, we also 
took into account the excitation of the modes in the frequency 
range [5,9] d"' . 

3. 1 . Fitting the dominant radial mode frequency 

In an initial step, we screened each evolutionary track 
(M, X, Z, Q-ov) from the ZAMS to the TAMS to find two consecu- 
tive models whose radial modes of a particular order encompass 
the observed frequency of 5.48694 d ' . We considered all pairs 
of consecutive models that explain the dominant mode as the 
fundamental or up to the fourth overtone. 

For each of those five paired models along the tracks that ful- 
fil this basic requirement, we subsequently computed additional 
models with a smaller time step to produce one model fitting 
exactly the frequency of the dominant radial mode, thereby fix- 
ing the age, or equivalently, the effective temperature, for that 
combination of (M,X,Z,ao^,n). In this way, we ended up with 
some 55000 models matching the dominant frequency. These are 
roughly equally distributed over all the considered radial orders 
n. 

For the seismically modelled /? Cep stars with a radial mode 
reported in the literature, this mode was either assumed to be the 
fundamental or first overtone, or higher overtones were ruled out 
because they were not compatible with the effective temperature 
or the gravity of the star Here, we are dealing with a star whose 
spectroscopic gravity was determined to be lower than the one 
of all those previously modelled stars (log g - 3.45). Hence, we 
did not place any limits on the overtone when interpreting of the 
dominant mode at this stage because we might be dealing with a 
higher overtone if the star is more evolved than usual for yS Cep 
stars. 

3.2. Requesting mode excitation for tlie dominant mode 

The oscillations of the ySCep stars are in general understood 
in terms of a heat mechanism active in the partial ionisation 
zones of iron-group elements (Moskalik & Dziembowksi 1992). 
While the majority of detected and well-identified low-order p 
and g modes is predicted to be excited by non-adiabatic os- 
cillation computations, some excitation problems occurred and 
remain to be solved for a few well-established observed non- 
radial mode frequencies, notably the ^ = 1, p2 modes of the stars 
vEri (Pamyatnykh et al. 2004) and yPeg (Handler et al. 2009, 
Zdravkov & Pamyatnykh 2009). To excite these modes in the 
models, one would need to increase the opacities of iron-group 
elements. In addition, the predicted frequency of this /f = 1, p2 
mode for appropriate models was found to be shifted with re- 
spect to the observed value for both vEri and yPeg. A similar 
excitation problem occurred for the star 12 Lac (Dziembowski 
& Pamyatnykh 2008; Desmet et al. 2009). 

Theory explains the excitation of the observed radial and pi 
non-radial modes well for the case studies of seismically mod- 
elled j0 Cep stars, except for the very massive (M ^ 24 M©) 09 V 
pulsator HD 46202, for which none of the detected frequencies 
is predicted to be excited (Briquet et al. 201 1). In particular, the 
detected radial modes of all seismically modelled /3 Cep stars of 
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Fig. 1. Models with predicted oscillation frequencies fitting the nine detected frequencies and fulfilling the mode identification as 
listed in Table[T] In the four upper panels, the dominant mode is the radial fundamental while in the lower left panel it is the first 
overtone and in the lower right panel the fourth overtone. Full coloured lines indicate excited modes while dashed coloured lines 
indicate stable modes. 



spectral type B agree with excitation predictions for these partic- 
ular modes. We are therefore justified in insisting that the dom- 
inant mode of HD 180642 be excited in present-day theoretical 
computations. To use this requirement as an additional constraint 
in our stellar modelling, we performed non-adiabatic computa- 
tions for the radial modes of all the ~55 000 selected models dis- 
cussed in the previous section with the code MAD (Dupret 2001, 
Dupret et al. 2002). 

For ~18 850 of the models that fit the frequency 
5.48694 d"' as a radial mode, this mode is predicted to be ex- 
cited. The distribution over the various overtones is as follows: 
9964 models have the dominant mode excited as fundamental, 
6681 as first overtone, 892 as second overtone, and 1324 as 
fourth overtone. There are no models with the third overtone ex- 
cited while there are with an excited fourth overtone. This may 
seem strange at first sight, but is explained by the occurrence of 
a resonance of the type Wp^ ^ 2tijp, . It was already known that a 
resonance ojp^ ^ Swp, occurs in polytropic models (Van Hoolst 
1996). Near-resonances were also found in that study, in partic- 
ular the one we also find here on the basis of more realistic mod- 
els (Fig. 9 in Van Hoolst 1996). In all the models with the fourth 
overtone excited, covering the mass interval [16.3, 20] M©, we 



also found the fundamental mode at 2.74 d" (and quite often 
also the first overtone) to be excited. This frequency does not oc- 
cur among the 127 dominant frequencies of HD 180642. These 
high-overtone models are more evolved and in closer agreement 
with the spectroscopic log g than the models with the fundamen- 
tal mode interpretation. 

3.3. Matching the additional non-radial mode frequencies 

In the subsequent steps of the modelling procedure, we could not 
rely on a unique {£, m) identification of the modes corresponding 
to the observed frequencies, as only limited additional informa- 
tion is available (see Table[Tli. Moreover, the additional modes 
all have much lower amplitudes than the one of the dominant 
mode, thus we did not prefer one above the other for the match- 
ing procedure. Thus we requested at once that all eight additional 
frequencies listed in Table[T|were fitted by model frequencies. 

We also imposed that the identification of the degree of 
the mode with frequency 7.35867d"' is fulfilled (^=0 or 3). It 
turned out that this frequency is not a higher overtone radial 
mode. For the frequency 8.40790 d"' , we requested it to have 
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one of the seven {{, m) combinations listed in Table[T] For each 
of those possibilities for {I, m), we used the rotational splitting 
along with the Veq value listed in Table[T]and model radius to de- 
duce the shift in mode frequency implied by rotation, according 
to the formula 



vnem = vnfo + \m /3„( X 86 400 d"' (1) 



(e.g., Aerts et al. 2010, Chapter 3), i.e., we limited the fre- 
quency matching procedure to first-order effects in the rota- 
tion frequency, which is consistent with our grid of spherically- 
symmetric evolutionary models. 

In the matching procedure, we screened all the models which 
excite and fit the dominant radial mode as described in the pre- 
vious section and fulfilling all the information in Table[T] for 
the eight additional measured frequencies, after shifting them 
according to Eq.([T]) and Table[T] The quality of the frequency 
matching was evaluated by the computation of a reduced 
statistic, where we used 0.02 d"' as the error for all the eight 
frequencies because it is a typical uncertainty on theoretically 
predicted frequencies for non-rotating main sequence pulsators 
due to differences in the input physics of the models (e.g., Moya 
et al. 2008). We then considered only those models with;i'^ < 1 
and for which each individual frequency is fitted to better than 
0.04 d ' . The latter constraint was adopted to avoid a rather 
large mismatch leading to < 1. 

In total, 2541 models survived our stringent matching con- 
straints, covering a mass range of [8.1, 19.9] Mq. Among those, 
some 1000 predict that the dominant frequency is the fundamen- 
tal mode, some 1200 that it is instead the first overtone and a few 
hundred that it is either p3 or ps . 

Six prototypical examples of schematic frequency spectra 
are shown in Fig.[T] The four that fit the dominant mode as 
the radial fundamental, predict that the 8.4 d ' frequency is an 
{ = 2,m = +1 mode and that the star has an equatorial rotation 
velocity of 26 km s ' . In the case of the first overtone model for 
the dominant mode, the 8.4d"' is an = +m - 1 mode for an 
equatorial rotation of 1 19 km s ' , and in case of the fourth over- 
tone model, the 8.4d ' is an ^ = +m = 2 mode for a rotation 
velocity of 52 km s"' . As can be seen from Fig.[T] all those cases 
give a good explanation of the observed frequencies. We stress 
that these are only six out of many good solutions. 

We point out that our procedure was not ideal for the three 
options for 8.40790 d"' with the highest equatorial rotation ve- 
locity listed in Table[T] for which second-order rotational effects 
may come into play at a level that can reach values above the the- 
oretical uncertainty we adopted. The estimate of the rotation pe- 
riod from the magnetic field study by Hubrig et al. (201 1) would 
have allowed us to drop these three options from the start. We 
preferred not to do so, however, as the derivation of the rota- 
tion period by Hubrig et al. relies on the assumption of dealing 
with a rigid rotator model and a centred dipole. While this is a 
very plausible assumption, and certainly the most obvious one 
to make, we wanted the seismic modelling to be as independent 
as possible of this result. This implies that our modelling is not 
optimal for those three options with high Veq, because we did not 
consider second-order effects in the rotation when computing the 
evolutionary models and the pulsation frequencies. Improving 
this would requke a full two-dimensional treatment of the equi- 
librium models and their frequencies, which is beyond the scope 
of this paper. 




Teff (K) 

Fig. 2. Kiel diagram showing all the models with a radial mode at 
the frequency 5. 48694 d"' excited and matching the additional 
eight higher mode frequencies listed in Table[T] The 1 cr error box 
for the spectroscopically derived T^ff and log g is shown as well. 
The dashed vertical lines indicate the restriction we imposed on 
the models. 
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Fig. 3. HR diagram showing the five models listed in the upper 
part of Table|2]as crosses on their evolutionary track (full lines). 
These models have excited oscillation modes whose frequencies 
match all the observed ones, except for the frequency 8.77 d , 
which fits but is predicted to be stable. The 38 dotted tracks each 
contain one model with excited oscillation frequencies match- 
ing the observed ones, except for the two modes with frequency 
above 8 d , which are predicted to be stable. 

3.4. Constraints imposed by spectroscopy 

There has been good agreement between spectroscopically de- 
rived effective temperatures and abundances, and their seismi- 
cally derived counterparts, from the method of frequency match- 
ing that we adopted here. This is not necessarily the case for 
the gravity of the seismically modelled y6 Cep stars. While good 
agreement was found for several stars (e.g., HD 129929, Dupret 
et al. 2004; 12 Lac, Desmet et al. 2009), the spectroscopic log g 
turned out to be 0.15 dex higher than the seismic log g for the 
star eOph (Briquet et al. 2007). For HD 46202, the discrepancy 
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even amounted to 0.24 dex (Briquet et al. 201 1). Given these pre- 
vious discrepancies, we prefer not to impose the spectroscopic 
gravity as a secure constraint in the modelling; this would imply 
that the dominant mode is a high overtone (n > 3) while we con- 
sider this unlikely. Similarly, the cuiTent metallicity of the star 
measured at its surface is not necessarily the initial one when the 
star was born, as effects such as atomic diffusion or unknown 
mixing may have occurred. Thus, we take a conservative atti- 
tude and request as an additional constraint from spectroscopy 
that the models must fulfil the Icr error bar in the effective tem- 
perature only. 

In Fig. 121 all 2541 models surviving the previous section are 
shown, along with the spectroscopically determined error box 
in the Kiel diagram. As can be seen, only two models would 
survive, if we were to impose the Icr error bar for the gravity 
of the star. For these two models, none of the observed non- 
radial modes is excited (see the following section). We consider 
this too restrictive and continued with all the models that fulfil 
Teff E [23 500, 25 500]K. There ai-e 357 models left after this 
requirement, with a mass coverage from 9.5 Mq to 15.6 Mq. This 
set of models still covers the entire considered range of values 
for the metalhcity, hydrogen fraction, and core overshooting. 

3.5. Excitation of tlie non-radial modes 

We subsequently checked the mode excitation of the low-order 
zonal frequencies for the 357 surviving models with the code 
MAD. Depending on the model parameters, we found excitation 
for the first radial overtone near frequency 7.1 d"' for masses 
between 9.7 and 13 Mq, for the £ - 1, pi mode with frequency 
in [6.0, 6.5] d"' for masses between 9 and IIM©, and for the 
^ = 2 f mode with frequency near 6.9 d ' for masses between 
9.4 and II.6M0. For f - 3 and f = 4-, the gi modes also have 
their frequency in [6.0,6.5] d ' and are found to be excited for 
masses between 9 and 12Mq while the ^ = 3 f mode has its 
frequency near 7.35 d"' and is excited for masses between 9.7 
and I2.2M0. These rough conclusions about the mode excitation 
are in qualitative agreement with the results on mode excitation 
found by Dziembowski & Pamyatnykh (2008), but these authors 
considered only a few representative models that have frequency 
spectra similar to those observed for the ySCep stars vEri and 
12 Lac, while we considered an entire grid of models. 

To constrain the parameters space for HD 180642, we first 
eliminated the models for which none of the non-radial modes 
with frequency above that of the dominant radial mode is ex- 
cited. This reduced the number of models to 221 and led to an 
upper limit in mass of 14.1 Mq. In particular, this removed the 
two models with the fourth overtone as a fit to the dominant fre- 
quency. Only five models able to predict that the dominant mode 
is the first overtone remained. All of these have only one non- 
radial mode excited in the frequency range [6,9] d ' . In each 
case, this corresponds to the ^ = 2, pi mode near 6.2 d ' . This 
is insufficient to explain the observed frequency spectrum and 
we are thus left with models whose fundamental radial mode fits 
the dominant frequency. 

Among the remaining 221 models, we found five that excite 
all but the highest frequency listed in Table[T] The properties of 
these five models are listed in the upper part of Table|2l These 
models are indicated with a cross in Fig. [3] and fulfil the spec- 
troscopically determined effective temperature. The frequency 
spectra of the two models whose mass is indicated in italics in 
Table|2]were compared with the observed ones in the middle 
panels of Fig.[T] Given that excitation problems occur for some 
modes with frequencies above 7d 'in the three well-known 



bright class members vEri, 12Lac, and yPeg, it is quite en- 
couraging to find models whose oscillation spectra are in full 
agreement with all the numerous observational constraints, ex- 
cept for the spectroscopic gravity and the excitation of the ^ = 3, 
pi mode. 

In addition to the five remaining models, there are 38 mod- 
els for which the mode excitation is fulfilled, except for the two 
highest frequencies in Tablefl] i.e., for one frequency less than 
the five models in the upper part of Table|2] The evolutionary 
tracks to which these also "acceptable" models in the lower part 
of Table|2]belong are indicated as thin dotted lines on Fig. [3] for 
reference with respect to the five best models. These 43 mod- 
els together cover a mass range of [10.5, 12.0] M©, the full range 
X € [0.68,0.74] of the grid, Z e [0.014,0.018], and a^^ be- 
tween zero and 0.2. That the initial internal metallicity of the 
star seems to be somewhat higher (Icr level) than the observed 
one at the surface for all the acceptable models, might have been 
introduced by the discrepancy in the gravity, which is not treated 
independently from the microtubulence in a spectroscopic analy- 
sis. It might also be due to effects of diffusion given HD 180642's 
longitudinal magnetic field of up to 700 G (Hubrig et al. 2009, 
2011). 

Almost all the seismic models have log g ^ 3.83, which is 
0.38 dex (2.5cr) above the best estimate from a classical spec- 
troscopic analysis (Paper II). We thus face another case where 
the seismic gravity differs from the spectroscopic gravity, just 
as for the ySCep pulsators eOph (B2IV, Briquet et al. 2007) and 
HD 46202 (09V, Briquet et al. 2011). Combining the seismic 
gravity with the equatorial surface rotation velocity of 26 km s"' 
and a radius ofR^ 6.8 R©, leads to a surface rotation frequency 
of 0.075 d (rotation period of 13.3 d). We thus find that a fac- 
tor of four occurs between the surface rotation frequency and 
the lowest frequency that was jointly detected in the CoRoT 
light curve and the spectroscopic time series. The frequency 
0.89870d"' listed in Table[T]is 12 times the rotation frequency, 
within the measurement errors. The seismic rotation period that 
we obtained from our best-fitting models is in excellent agree- 
ment with the most probable period deduced from the magnetic 
field measurements. 

3.6. The full oscillation spectrum 

Continuing further with the models in Tabled we considered 
their entire oscillation spectrum in the range [0.2,40] d"' , tak- 
ing into account all rotational splittings according to the Ledoux 
formulation (e.g., Aerts et al. 2010) for a rotation frequency 
of 0.075 d"' , i.e., assuming rigid rotation. In Fig.lH we show 
the mode inertia (defined as Eq. (3.140) in Aerts et al. 2010) 
versus the oscillation frequencies for the seismic model with 
M = II.4M0 Hsted in Table|2] in the range [0,20]d ', for 
modes with i = 0, . . . , 4. We note that the inertia follow a 
smooth curve as the radial order increases, except for some spe- 
cific modes of ^ = 3 and 4 where the details of the trapping 
in the interior change owing to the interplay between the mode 
frequency and the Briint-Vaisala and Lamb frequencies. The de- 
tected frequencies in the best fitting CoRoT light curve model 
deduced in Paper I are shown as dotted vertical lines, except for 
the harmonics of the dominant frequency. Full symbols denote 
modes predicted to be excited from the non-adiabatic computa- 
tions done with the code MAD, while the open symbols are pre- 
dicted to be stable. 

All the combination frequencies (listed in Table 2 of Paper I) 
in the range [4, 8.5] d ' can be closely identified with frequen- 
cies of excited modes. This is also true for the four frequencies 
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Table 2. Stellar models in agreement with the observational properties of HD 180642. 



Mass 


Z 




X 




log S 


log(L/Lo) 


R 


Age 


X, 


(Mo) 








(K) 


(cgs) 




(Ro) 


(lO^r) 




11.4 


0.018 


0.00 


0.68 


24520 


3.829 


4.177 


6.80 


0.1243 


0.209 


11.6 


0.018 


0.05 


0.68 


24880 


3.832 


4.206 


6.84 


0.1220 


0.233 


11.7 


0.018 


0.05 


0.68 


25000 


3.834 


4.217 


6.86 


0.1200 


0.235 


11.7 


0.018 


0.05 


0.70 


24550 


3.834 


4.186 


6.86 


0.1296 


0.246 


11.8 


0.018 


0.05 


0.70 


24670 


3.835 


4.197 


6.88 


0.1275 


0.248 


10.5 


0.014 


0.05 


0.68 


23880 


3.814 


4.111 


6.64 


0.1482 


0.182 


12.0 


0.018 


0.05 


0.72 


24460 


3.837 


4.187 


6.91 


0.1332 


0.262 




Fig. 4. Mode inertia for the oscillations of the model with M - 11 .4 M©, whose other parameters are listed in Table|2l Thick crosses 
indicate the radial fundamental and first overtone, which are both excited. Thin crosses indicate the higher overtone radial modes, 
which are not excited. Filled symbols indicate excited non-radial modes and empty symbols stable modes according to the indicated 
label. The meaning of the symbols is as follows: cross: ( - Q, squares: € - \, circles: { - 2, downward triangle: f - 3, upward 
triangle: { - 4. For clarity, we show only the excited g modes in the region below 1.3 d"' in the large panel, while all g modes are 
shown in the inset. The dotted vertical lines indicate the frequencies detected in the CoRoT light curve (Paper I), nine of which have 
also been detected in ground-based data (Paper II). 



in the range [0.8, 1.2] d ' , which can be explained by excited 
g modes of intermediate to high radial order (n typically be- 
tween 10 and 30). All other combination frequencies found in 
Paper I can be attributed to model frequencies of non-excited p 
modes taking into account the observational and theoretical un- 
certainties. If the frequency 0.3 d ' were produced by pulsation, 
it would have to correspond to an = 3 mode (Paper II). There 
are indeed various £ = 3 modes available, taking into account 
the rotational splitting (see the inset of Fig.|4]l. These have radial 
orders from 53 to 60. Given the magnetic field detection for the 
star, it might also be that this frequency detected in both the pho- 
tometry and spectroscopy originates from low-amplitude surface 



inhomogeneities. We are unable to discriminate between these 
two options. In any case, we interpret the detection of all these 
modes with combination frequencies that are not excited by the 
heat mechanism in terms on a non-linear resonant mode excita- 
tion by the dominant radial mode, as already proposed in Paper I. 
That this mechanism seems to be at work and results in frequen- 
cies beyond the usual frequency regime of the /3 Cep stars is sup- 
ported by the low inertia of these modes (Fig.HJl, which makes 
it much easier to excite them than low-order g modes. The high 
amplitude of the radial mode thus seems to trigger the excitation 
of p modes at frequencies between [8.5, 20] d"' , which would 
otherwise be stable. 
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Fig. 5. Echelle diagrams of the CoRoT frequencies for 
HD 180642 for the large spacing determined by Belkacem et 
al. (2009, upper panel) and in Paper I (middle panel). For com- 
parison, we also show the echelle diagram of the 08.5V star 
HD 46149 reproduced from Degroote et al. (2010b), where clear 
ridges produced by stochastically excited modes are found. The 
size of the symbols scales with the power of the frequency. 




1.5 2 
Frequency difference (d ^) 

Fig. 6. Fraction of the occurrence of values for the frequency 
differences between zonal p modes of degree { - 0, . . . , 4 in 
the frequency regime [1 1, 26] d"' compared with the spacing of 

I. 17 d ' reported by Belkacem et al. (2009). The full line is for 
the model of 1 1 .4 M©, while the dotted line is for the model with 

II. 7 Mq indicated in italics in Table|2] 



The results discussed for the model with M - 1 1 .4 re- 
main qualitatively the same for the four other models listed in 
the upper part of Table|2] 

3. 7. Revisiting fhe case of solar-like oscillations 

We used the set of models that fulfil all the seismic require- 
ments of the p Cep-type modes, as well as the combination fre- 
quencies, to revisit the claim of solar-like oscillations in the star 
(Belkacem et al. 2009). In Paper I, the authors already pointed 
out that the value found for a possible frequency spacing de- 
pends on the prewhitening procedure adopted to model the light 
curve. Moreover, it was found in Paper! that several of the 
modes with combination frequencies considered by Belkacem et 
al. (2009), who made their computations for the frequency range 
[1 1, 26] d"' , are phase-locked, which is not what is expected for 
stochastic forcing. The interpretation of the residual light curve 
in terms of solar-like oscillations was also questioned by Balona 
et al. (2011), in view of the lack of such oscillations in Kepler 
data of B-type pulsators, because these data have higher preci- 
sion than the CoRoT data. One expects solar-like oscillations 
in stars with similar fundamental parameters than HD 180642 
to have comparable amplitudes, and these amplitudes would be 
far above the detection limit of Kepler data. In the sample of 
Balona et al. (2011), seven stars have an effective temperature 
above 20 000 K and are positioned in the joint part of the yS Cep 
and SPB instability strips, but only one of these seven has a tem- 
perature close to the one of HD 180642. This star's variability 
is probably due to rotation and at the mmag level, i.e., a factor 
of more than ten below the one of our target, so the comparison 
between the Kepler stars and HD 180642 in terms of solar-like 
oscillations is premature. 

On the observational side, we made an extensive search for 
ridges in echelle diagrams constructed for the CoRoT resid- 
ual light curve of HD 180642 discussed in Paper 1, adopting the 
same method as in the case of the CoRoT data of the 08.5V 
pulsator HD 46149 (Degroote et al. 2010b). In doing so, we con- 
sidered numerous values for a large frequency spacing. We failed 
to identify any ridges, in particular also for the value of the spac- 
ing reported in Belkacem et al. (2009). This model-independent 
result stands in contrast to the clear ridges present in the echelle 
diagram of HD 46149. A comparison between the echelle dia- 
gram of HD46149 and two for HD 180642, one based on the 
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large spacing by Belkacem et al. (2009) and the most likely one 
reported in Paper 1, is shown in Fig.|5] The recently discovered 
solar-like oscillations in the frequency regime above that of the 
heat-driven modes in a (5 Set star observed with the Kepler satel- 
lite did lead to clear ridges in the echelle diagram of that star 
(Antochi et al. 2011), as is the case for HD 46149. 

We then considered our seismic models that explain all the 
modes of HD 180642 in the frequency regime [0, 20]d 'in 
terms of heat-driven modes with non-linear mode couplings 
(listed in Table|2]i. Irrespective of the cause of the power in the 
range [11,26] d"' , we first computed all frequency differences 
between the axisymmetric mode frequencies of consecutive ra- 
dial order for the same degree £ in that interval and compared 
all these spacing values with the one reported by Belkacem et 
al. (2009), with an allowed uncertainty of 0.02 d ' . For the 
model with M - II.4M0, this spacing occurs only once, i.e., 
between the £ = 3,m = 0, gi and g2 modes. It also occurs 
once between the ^ = 3,ot = 0, p3 and p4 modes of the model 
with M - II.8M0, while it never occurs for the three other 
models in the upper part of Table|2l Subsequently, we com- 
puted all frequency differences for all the axisymmetric modes 
of ^ = 0, ...,4 in the frequency interval [ll,26]d"' , to test 
whether any particular spacing occurred. The outcome of these 
computations is shown in Fig.|6]for the two models indicated in 
italics in Table|2] It can be seen that no frequency spacing can be 
discerned among the model frequencies. 

4. Concluding remarks 

We have managed to explain the extensive observational seismic 
input deduced from CoRoT and ground-based photometry along 
with high-resolution spectroscopy of the /3Cep star HD 180642 
by means of present-day models. All of the 30 detected frequen- 
cies are fitted by model frequencies, of which 16 are predicted to 
be excited and 13 are combination frequencies that can be linked 
to the dominant radial fundamental mode with an extremely 
large amplitude for this type of pulsator. Our modelling resulted 
in a mass between 11.4 and 11.8 Mq and a zero or very low core 
overshooting parameter, along with an age of between 12 and 13 
million years (corresponding to a central hydrogen fraction be- 
tween 0.21 and 0.25) if we allow an excitation problem for only 
one of the independent oscillation modes, namely for the £ - 3, 
pi mode with a frequency 8.77 d . This is a frequency regime 
where similar excitation problems have been found for this type 
of pulsator (e.g., Dziembowski & Pamyatnykh 2008, Zdravkov 
& Pamyatnykh 2009). All other observed frequencies that are 
not excited by the heat mechanism can be explained as eigen- 
modes triggered by the large-amplitude radial mode by means 
of non-linear resonant mode excitation. 

We are unable to support the findings of Belkacem et al. 
(2009) from the seismic models which explain the oscillation 
spectrum ofHD 180642 in the frequency regime [0,20] d ' , i.e., 
the nature of the power excess in the p-mode frequency regime 
of HD 180642 is different from the one of "classical" solar-like 
oscillations caused by the stochastic forcing due to convective 
motions. We instead attribute the occurrence of time-dependent 
amplitudes and phases in the frequency regime above the radial 
fundamental mode to the effect of the non-linear resonant mode 
coupling which is activated by the large amplitude of the dom- 
inant mode. This provides a natural explanation for the lack of 
detections of solar-like oscillations in B stars monitored by both 
the CoRoT and Kepler missions. All of the B pulsators observed 
by these satellites indeed have much lower amplitudes (at least 
by a factor ten below the dominant mode of HD 180642) so reso- 



nant mode coupling is not expected for these, while solar-Uke os- 
cillations would have been detected if they existed in such types 
of stars. 

We find, once more, a discrepancy between the spectroscopic 
log g and the seismic log g, as was found before in two other 
seismically modelled OB stars. We assign this to the system- 
atic uncertainty when deducing log g from the wings of spectral 
lines of one snapshot spectrum or from an average of a time 
series of spectra. In both cases, the pulsational broadening is 
not taken into account in an appropriate way. Indeed, oscilla- 
tions change the shape of spectral lines in an asymmetric and 
time-dependent way, i.e., the level of deformation depends on 
the phase in the overall oscillation cycle. To compensate for this 
deformation and absence of pulsational broadening, one usually 
introduces a time-independent ad-hoc parameter called macro- 
turbulence (e.g., Aerts et al. 2009) and/or a rotational profile 
with a wrong vsin /, each of which are usually deduced from 
metal lines and have a different functional shape than the true 
time-dependent line-broadening function. This is then compen- 
sated for by an incorrect value of the gravity and/or microtur- 
bulence, which can either be too high or too low. In addition 
to this, the normalisation of the spectra, which is a necessary 
manipulation of the data required to deduce the position and 
shape of the line wings with respect to the continuum, is by 
no means a trivial operation and in itself introduces systematic 
uncertainty. In the present case of HD 180642, for instance, the 
gravity was deduced from four Balmer line wings assuming a 
rotational profile of 44 km s"' , while a problem was encountered 
in deriving the microturbulence in a consistent way, because a 
discrepancy of 5kms"' occurred between the values deduced 
from the N and O lines for this quantity. Our seismic analysis in- 
ferred that V sin / - 25 km s ' only. There are thus various uncer- 
tainties to be considered for quantities that are not independent. 
The seismic gravity is a quantity determined by the entire stellar 
structure, which is very strongly constrained for HD 180642 by 
the numerous seismic properties measured for this star. In view 
of the discrepancy, we advise the use spectroscopically derived 
Zo^ ^-values of pulsators with a 3cr error bar Temperature es- 
timates are not affected by this uncertainty as they rely on an 
ionisation balance of Si for which three stages are available in 
the case of /JCep stars, hence this quantity is not very dependent 
on the line shapes. 
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